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I. INTRODUCTION 

Describing the relation between particle and energy 
currents is at the heart of thermoelectrics.^^® For dc driv¬ 
ing with small temperature gradients and bias voltages, 
linear-response relations between the currents and the 
applied forces constitute the basis to describe thermo¬ 
electric phenomena. When combined with the principles 
of thermodynamics, the resulting theory has the beauty 
of simplicity and the strength of high predictive power. 
Specifically, it allows for a successful characterization of 
the efficiency of various thermoelectric machines in terms 
of the figure of merit introduced by Ioffe in 1949.® 

An important challenge is to incorporate genuine quan¬ 
tum effects associated with coherent transport in nano¬ 
devices into this theoretical framework for thermoelec¬ 
tric effects. Here, we address how to include adiabatic 
quantum pumping as a paradigm of coherent-transport 
effects into a suitably generalized thermoelectric frame¬ 
work and explore the fundamental relations of the cor¬ 
responding quantum machines. Quantum pumping gen¬ 
erates nonzero dc currents by locally applying purely ac 
drivings to a quantum coherent conductor. It gen¬ 
erates both charge and energy currents^^, enables heat 
pumping, and the exchange of work between different 
driving forces^^’^®. The aim of the present work is to ex¬ 
tend the linear-response theory of thermoelectric effects 
to systems under adiabatic driving. To this end, we need 
to include the energy flux between the electrons and the 
ac forces on an equal footing with the heat and particle 
fluxes. 

Figure 1 shows the setup that we have in mind. It 
consists of a central coherent conductor which is cou¬ 
pled to two reservoirs. In conventional thermoelectrics, ^ 
the two reservoirs differ in both temperature and chem¬ 
ical potential. A thermal engine converts a temperature 
difference into electric power. As a consequence of the 
Second Law, the efficiency of this conversion process is 
limited by the Carnot efficiency rjc = {T 2 — Ti)/T 2 where 
T 2 > Ti denote the temperatures of the reservoirs. The 


optimal efficiency that can be reached for a specific device 
is controlled by its figure of merit or ZT value, ^ 


0 = 00 


y/l + ZT- 1 
^/1 + ZT+l' 


( 1 ) 


The ZT value can be expressed in terms of the linear- 
response coefficients of the device, relating charge and 
heat currents to bias voltage and temperature gradient.^ 
The thermal engine can also be operated in reverse, 
realizing a refrigerator which invests electric power to 



FIG. 1. Sketch of the setup. A coherent quantum conductor 
is driven by time-periodic potentials and connected to two 
reservoirs biased by (a) a chemical-potential difference Sp, or 

(b) a temperature gradient ST, or both. Charge A_r, heat 

and power W are exchanged between the reservoirs and the 
ac sources. The solid (dashed) arrow indicate (a) the motor 
(generator) mode of the device, and (b) the heat engine (heat 
pump) mode. 

























2 


continuously extract heat from the colder reservoir. The 
maximal efficiency of this device is given by the appropri¬ 
ate Carnot efficiency ijc = Ti/{T 2 — Ti) and in terms of 
this Carnot efficiency, the optimal efficiency for a specihc 
device is again given by Eq. 1.^ 

In this paper, we consider setups in which the coherent 
conductor is subject to a set of ac potentials in addition. 
For definiteness, we will consider reservoirs which have ei¬ 
ther different chemical potentials [Fig. 1(a)] or different 
temperatures [Fig. 1(b)], although our theory could read¬ 
ily be applied to situations which combine ac potentials 
with both chemical-potential and temperature gradients. 
The physics of these setups can be understood by analogy 
to the Archimedes device, a pipe with a rotating screw, 
which can be used to pump water against gravity. This is 
a classical analog of an adiabatic quantum pump, where 
ac driving pumps a certain amount of electric charge per 
cycle. Specifically, this charge can be pumped against an 
applied dc bias voltage,in which case quantum pump 
realizes a generator. 

The Archimedes screw can also be operated in reverse, 
with water flowing between the reservoirs by gravity and 
setting the screw into rotational motion. An analogous 
effect can be used to turn an adiabatic quantum pump 
into an adiabatic quantum motor. This is most easily 
understood when imagining that the time dependence of 
the ac potentials derives from the (classical) dynamics of, 
say, one or more mechanical degrees of freedom. Then, 
a charge current pushed through the coherent conductor 
will set the mechanical degrees of freedom into motion. 

Generator and quantum motor are driven by a bias 
voltage and correspond to the setup sketched in Fig. 1(a). 
Alternatively, we can also consider devices involving tem¬ 
perature gradients instead of bias voltages, which realize 
heat pumps and heat engines. Such a device is depicted in 
Fig. 1(b). The devices in Fig. 1 are examples of nanomo¬ 
tors and nanoengines, which have received much atten¬ 
tion recently. We note that the effect of ac potentials 
on the conventional thermoelectric effects has been stud¬ 
ied in a number of recent papers. 

The aim of the present work is to extend the linear- 
response theory of thermoelectrics to such nanomotors 
and nanoengines, to understand their efficiencies, and to 
identify appropriate hgures of merit. This program poses 
several conceptual questions: (i) We need to identify the 
current that complements the charge and heat currents 
and accounts for the effects of the ac potentials. Simi¬ 
larly, we need to identify the affinity that complements 
the (scaled) temperature difference and bias voltage, (ii) 
We need to develop the generalized linear response theory 
which includes these additional quantities. While this is 
a conventional linear-response theory for traditional ther¬ 
moelectrics, the ac potentials are not actually weak but 
only slowly varying, (iii) We finally need to identify ap¬ 
propriate efficiencies and figures of merit. We will see 
that the latter also differ in essential ways from those 
defined in conventional thermoelectrics. 

In Sec. II, we generalize linear-response theory to in¬ 


clude the response to the adiabatically varying ac poten¬ 
tials in addition to the applied bias voltage. We do this 
by working to linear order in the rate of change (or ve¬ 
locity) of the ac potentials. We find that this can be 
done in a manner which closely resembles the deriva¬ 
tion of Kubo formulas in linear response theory. Conse¬ 
quently, we derive general Kubo-like expressions for the 
response of both the charge current and the generalized 
forces conjugate to the ac potentials. These expressions 
imply that the response coefficients satisfy Onsager-like 
relations and are thus not independent of one another. 
In Sec. Ill, we generalize the thermodynamical frame¬ 
work to include the time-averaged work per unit time 
performed by the ac forces as a third flux, along with 
the heat and particle fluxes. We also identify the scaled 
frequency huj/T of the driving as the appropriate third 
affinity, complementing the temperature and chemical- 
potential differences. In Sec. IV, we dehne and analyze 
efficiency and hgure of merit for the various quantum 
machines sketched in Fig. 1. We find that the dehni- 
tion of the appropriate figure of merit analogous to the 
ZT value differs in characteristic ways, reflecting the fact 
that the usual off-diagonal thermoelectric response coef¬ 
ficients, the off-diagonal coefficients involving the third 
flux or affinity do not enter into the entropy production. 
To illustrate these concepts, we apply our theory to an ex¬ 
ample device in Sec. V which is based on a simple model 
for a quantum pump. We summarize in Sec. VI. 


II. ADIABATIC RESPONSE AND ONSAGER 
RELATIONS 

We begin by evaluating the forces and currents induced 
by a set of time-periodic parameters in the adiabatic ap¬ 
proximation. We will see that this can be done in close 
analogy to linear-response theory, allowing us to derive 
Onsager-like relations. 

We collect the parameters Vi{t) of the Hamiltonian % 
into a vector V(t) = V(t-|-T) = V 2 {t ),...) so that 

a = 'H{'V{t)), where T = 27r/a; is the driving period. 
Quite generally the Hamiltonian of the system can be 
expressed as 


H(V(t))=?io-E^:>^^W> (2) 

3 

where Hq is the time-independent part of the Hamilto¬ 
nian and Fj are hermitian operators that play the role of 
generalized forces 


F(t) = 


diiit) 
dY{t) ■ 


( 3 ) 


The quantum expectation value tr{pF} in terms of the 
electronic density matrix p are just conventional forces 
when the Vj denote regular cartesian coordinates of a 
classical system obeying Newtonian dynamics. 
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At lowest order in the adiabatic approximation, the 
system is described by the frozen density matrix p* for 
the Hamiltonian "Ht with t treated as a parameter. Ac¬ 
counting for the temporal variation of V (t) to lowest or¬ 
der, we can approximate the time evolution operator as 


U{t,to)-Texp{-iHt{t-to)-i [ dt'{t-t')F -Vit)}. 

Jto 

(4) 

To linear order in the small ’’velocity” V(t), we can now 
follow the usual steps of linear response theory^® and 
express the expectation value 0{t) of an observable O at 
time t as 


Jto 

= (d)* + Af^-V(t). (5) 


Here, the operators 0{t) and F(f') are defined in the 
Heisenberg representation with respect to the frozen 
Hamiltonian 'Ht and (.. .)t denotes the expectation value 
with respect to the frozen density matrix pt. The re¬ 
sponse function can be expressed through the re¬ 
tarded adiabatic susceptibility — t') = —i9{t — 

f')([0(t),F(t')])t. We now expand the frozen average to 
linear order in an applied bias yielding {0)t — 
where the linear-response coefficient Ap'’ is given by the 
usual Kubo formula. Applying this procedure specifically 
to the charge current J^(t) and the forces F(t) (and post¬ 
poning the heat currents and temperature gradients for 
further below), we obtain 


(j%t) U ^ 


( 6 ) 


to linear order in 6p and V(t). 

The terms in Eq. (6) have clear physical interpreta¬ 
tions. The first term on the right hand side collects 
the currents and forces evaluated with the frozen den¬ 
sity matrix pt in equilibrium (be., for 5p = 0). These 
terms have zero mean when averaged over one period 
of the ac fields. The forces can be thought of as con¬ 
servative Born-Oppenheimer forces and expressed as a 
gradient of the equilibrium energy of the system with 
respect to V(t). For several potentials this term may 
lead to exchange of work between the different forces Fj 
without dissipation. Such processes were considered in 
Refs. 12 and 13. Adiabatic quantum pumping of charge 
by the ac potentials is described by Aj-^, while A/^ cap¬ 
tures the modification of the forces by the applied bias 
6p. Both contributions are generally nonzero when av¬ 
eraged over a period, implying that this contribution to 
the force is raonconservative. This was discussed for non¬ 
interacting electrons coupled to adiabatic nanomechani¬ 
cal systems^®’^^ and nanomagnets.In the latter case, 
this corresponds to a spin-transfer torque. The diago¬ 
nal components describe the usual conductivity through 
A“ and the velocity-dependent force through A{^ . In 


time-reversal symmetric systems, the latter is symmetric 
and describes a frictional force. Without time-reversal 
symmetry, A{^ may have an antisymmetric part which 
is analogous to the Lorenz force. 

The derivation of the response coefficients Af follows 
standard linear-reponse theory, including the “adiabatic” 
response to the ac potentials. Consequently, it is natural 
to expect that the response functions Af satisfy Onsager- 
like relations. In fact, these can be derived in the usual 
manner, as shown in detail in App. A. Thus, we find the 
generalized Onsager relations 


AT{B)=AT{-B) , A"(R) = s.s,Aj/(-R) 

Af{B) = s,Af (-B), (7) 

where the sign Sj = ± depends on the parity of the op¬ 
erators Fj under time reversal. As the derivation of On¬ 
sager relations is very general, these relations are valid at 
finite temperature T and in the presence of many-body 
interactions.^® 

The second line in Eq. (7) imposes a relation between 
the adiabatic quantum pumping of charge (as described 
by A^^) and the nonconservative force (as described by 
AJ'’). This relation which is valid in the adiabatic regime 
was previously found for noninteracting adiabatic quan¬ 
tum motors at zero temperature and B = 0.^^ It has 
been pointed out that time-reversal symmetry - by way 
of Onsager-like arguments - does not imply symmetry of 
the pumped charge under magnetic-field reversal unless 
the system has additional spatial symmetries.The 
relation in Eq. (7) implies that there is still an Onsager 
relation associated with the pumped charge, but it does 
not relate the pumped charge to itself but rather to the 
nonconservative force in response to an applied bias. 

Before closing this section, we comment on how to in¬ 
clude heat currents and thermal gradients into this linear- 
response scheme. Within linear response, we can readily 
extend Eq. (6) into a 3 x 3 matrix equation 

6p 

ST 
V{t) 

( 8 ) 

Here, we can identify the thermal conductance A^® re¬ 
lating the heat current J^{t) to ST as well as the 
usual thermoelectric coefficients A‘^‘^ and A'^°. In addi¬ 
tion, our scheme includes the coefficients A^^ and A{'^ 
which describe the generation of heat currents by a time- 
dependent driving (quantum pumping of heat) and the 
generation of a nonconservative force in response to a 
temperature gradient, respectively. 

The treatment of a temperature gradient within the 
Kubo approach is less straightforward, but has been ad¬ 
dressed numerous times in the literature.An alter¬ 
native route is to calculate the relevant observables with 
a non-equilibrium technique, such as the Keldysh or scat¬ 
tering matrix formalisms, and to perform the expansions 
in ST, Sp, and V a posteriori. (This is the route which we 



/ nt) \ ( Jt\ ( Ar A“« Af \ 

= J« + Af Af Af 

\ m ; \ A^ Af^ A{f y 
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follow in Sec. V). Either approach yields the additional 
Onsager relations 

Kr{B)=Kl\-B) , Kl\B)=Kr{-B) 

Af{B) = s,^\-B) (9) 

complementing Eq. (7). The first line corresponds to the 
usual thermoelectric Onsager relations. The second line 
contains the additional Onsager relations relating pump¬ 
ing of heat current and the force generated in response 
to an applied thermal gradient. 

III. GENERALIZED THERMOELECTRIC 
FRAMEWORK 

Conventional thermoelectrics considers particle and 
heat currents in response to chemical-potential and tem¬ 
perature differences. In the presence of ac driving as 
in the devices in Fig. 1, we have to take into account 
the pumping of particles and heat as well as the work 
performed by or on the ac potentials on the same foot¬ 
ing. To develop the corresponding generalized thermo¬ 
electrics, we first consider the entropy production of the 
system. After averaging over one period of the ac driv¬ 
ing, the net dissipation occurs only in the electrodes and 
we can write 

+ (10) 

where the average heat flux in lead a is given by 

Qa = Ea - (J-aNa. (11) 

The energies and particle numbers Na satisfy the 
conservation laws 

Nii = -NL, EL+En = W. ( 12 ) 

While particle-number conservation takes the same form 
as in standard thermoelectrics, energy conservation must 
account for the additional work W performed by the ac 
potentials on the electron system. The corresponding 
power can be expressed as IE = — J2j yielding 

the entropy production 

S = + (13) 

3 

to linear order in the applied bias Sfi = fj,L — ^J-R and tem¬ 
perature difference 6T = Tl — Tn- Note that after aver¬ 
aging over a period, the conservative Born-Oppenheimer 
forces in Eq. (6) do not contribute to entropy production. 
Then, the power can be expressed in linear response and 
for (5r = 0 as 




Here, the first term on the right-hand side describes the 
work performed by the nonconservative force originating 
from the applied voltage (5/i {5T would contribute a sim¬ 
ilar term) and the second term is the dissipated power 
due to a frictional force on the ac potentials. 

In conventional thermoelectrics, one defines the par¬ 
ticle and heat fluxes Ji = Nr and J 2 = Qr as well 
as the corresponding affinities Xi = 5^/T and X 2 = 
^ 2^/212 41,42 -pQ extend thermoelectrics to the present sit¬ 
uation, we need to identify appropriate fluxes and affini¬ 
ties for the ac driving terms. 

At first sight, Eq. (13) may suggest to define the —Fj as 
fluxes and the Vj/T as the associated affinities. However, 
Eq. (13) holds only after averaging over one period. Be¬ 
fore time averaging, the conservation laws involve addi¬ 
tional terms^® and the forces Fj (t) contain contributions 
that are conservative. We can identify an appropriate 
affinity by noting that after averaging, the first term in 
Eq. (14) is proportional to w, while the second term is 
proportional to It is thus natural to define the affin¬ 
ity A 3 = huj/T with associated flux J 3 = lT/(/iaj).^^d 2 
Thus, Eq. (13) yields 

= ( 15 ) 

3 

for the rate of entropy production. 

We complete our quantum thermoelectrics scheme by 
linear-response relations between fluxes and affinities, 

J^ = Y,L^kXk. ( 16 ) 

k 

The linear-response coefficients Ly are readily related to 
the coefficients which appeared in Eq. ( 8 ). Indeed, we 
have 

Ln=TAf, Li2=r"Af, L13 = TA^^ • v, 

L21 = TAf, L22 = T^Af, L23 = TAf • V, ( 17 ) 
L31 = -TA{^ ■ V, L32 = -T^A{^ ■ V, 

L33 = -Tv^ • A*" • V, 

where we defined v through V = hujv and denotes 

the transpose of v. 

Thus, the coefficients Lij also obey Onsager relations, 
namely 

U,{B) = Lr{-B) , U,{B) = ±Lj,{-B), (18) 

with i ^ j. The sign in the second relation depends on 
the behavior of the fluxes under time reversal. Assuming 
time reversal from now on (and thus B = 0), this yields 
the relation L 12 = ^ 21 , which is well known from the 
usual theory of thermoelectrics, as well as L 13 = —T 31 
and L 23 = —L^ 2 - It is important to note that the off- 
diagonal response coefficients have the same sign in con¬ 
ventional thermoelectrics, while they have opposite signs 
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when either J 3 or is involved. We will see that be¬ 
low this has significant consequences for the definition of 
figures of merit for the devices in Fig. 1. 

The transport coefficients can be directly calcu¬ 
lated from the coefficients A, which are in turn given in 
terms of the susceptibilities Xt(a;). Another possibility is 
to start from the expressions for the charge, heat, and 
work currents, to perform the expansions in /ku, i5/r, and 
5T, and to identify the coefficients L from the resulting 
expressions. For noninteracting systems, this procedure 
is rather straightforward. 

In Sec. V, we will illustrate our general theory for a 
general noninteracting model of a two-terminal conduc¬ 
tor and evaluate the various response coefficients explic¬ 
itly. This will rely on Green function^^’^^do scatter¬ 
ing matrix^^d 2 expressions for the response coefficients 
which we derive by the procedure described in the pre¬ 
vious paragraph. Details of the calculations are given in 
App. C. The calculations start with the expressions for 
charge current [Eq. (B2)], heat current [Eq. (B4)], and 
work current [Eq. (B 6 )] for this model. Performing the 
expansions in hu, (5/i, and 5T, we find the explicit formu¬ 
las for the Lij given in App. D. One can also check that 
these expressions for the response coefficients satisfy the 
generalized Onsager relations Eq. (18), as they should. 


We are now ready to characterize the performance 
of adiabatically-driven quantum motors or generators in 
terms of efficiencies and figures of merit. The efficiency 
^mot gj: ^ motor is measured by the ratio of the work per 

unit time —W performed on the ac potentials and the 

power NnSfi/e injected by the voltage source. Similarly, 
the efficiency of the generator ry®®" is given by the inverse 
of this ratio, so that 


V 


mot 


1 


Tygen 


-W 

NnSfi/e 


( 20 ) 


Note that we have defined = ^J-R + (as well as 
Tl=Tr = T), cf. Fig. 1. 

We first show that the second law of thermodynamics 
implies an upper limit for these efficiencies. Using Eqs. 
( 10 ), ( 11 ), and ( 12 ), we find 


W = TS-^Nh- 

e 

Substituting this into Eq. (20), we obtain 


^mot ^ ^ 1 

NnS^l/e 


( 21 ) 


( 22 ) 


IV. EFFICIENCY AND FIGURE OF MERIT OF 
QUANTUM MACHINES 

A. Motors and generators 

Consider a situation with applied ac driving forces and 
a dc bias d/i, but uniform temperature T. The device in 
Fig. 1(a) can operate as a quantum motor or generator. 
When the ac potentials pump particles into the reservoir 
with lower chemical potential, the gain in electrical en¬ 
ergy can be used to perform work on the source of the ac 
potentials. This occurs for L^iS^/T < 0 and corresponds 
to a motor as the work performed on the ac potentials 
can be further transformed, say, into mechanical work.^^ 
When reversing the sign of Sfi and thus LsiSfi/T > 0, 
the ac potentials pump particles into the reservoir with 
higher chemical potential and we have a generator. 

Using X 2 = 0, the rate of entropy production becomes 


Now, the second law of thermodynamics demands S > 0. 
Moreover, the current flows with the potential drop 6fj, 

in the motor, Nji6^ > 0 , but against the potential drop 
for a generator, N< 0. Consequently, we find that 
both and 77 ®®" are upper bounded by unity. 

The efficiency in Eq. (20) can also be written as 


mot ^ 1 ^ _ ^ 3^3 

T^Sen XiJi 


(23) 


and the currents expressed through their linear-response 
expressions (16). Still assuming time-reversal symmetry, 
so that Li 3 = —T 31 , we can then maximize the efficiency 
as a function of Xi at fixed A 3 . One finds that the 
efficiency is maximized for 


L11L33 ± VTiiTsadetL 

Xi = - - — - -A 3 . (24j 

^llJ^13 


S = LnX^ + L33A32 + (Li3 + L31) ^1^3- ( 19 ) 

Interestingly, the last term on the right-hand side van¬ 
ishes due to the Onsager symmetry L 13 = —T 31 , and 
the coefficients L 13 and L 31 do not affect the entropy 
production. As a consequence, the second law of ther¬ 
modynamics imposes Ln > 0 and L 33 > 0. This is in 
contrast to conventional thermoelectric, where the off- 
diagonal response coefficients are symmetric, L 12 = T 21 , 
and do contribute to entropy production. In the latter 
case, the second law imposes detL = T 11 L 22 — L 12 > 0 
in addition. 


with +(—) for motors (generators). Alternatively, we can 
fix Xi is fixed and maximize the efficiency as a function 
of A 3 . In this case, one finds 

± y/LiiL^^detL 

A 3 = -----Ai, (25) 

T 33 T 13 

again with -l-(—) for motors (generators). Substituting 
Eqs. (24) and (25) into Eq. (23), we find for the maximal 
efficiency 


yr+c -1 

yr+c+i 


( 26 ) 
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and identify the figure of merit as 


ficiencies. Using Eqs. (10), (11), and (12), we find 


L 11 L 33 


(27) 


w = TrS + 


Tl - Tr 
Tl 


Ql- 


(29) 


Note that and C, are valid for both motors and gen¬ 
erators. 

Equations (26) and (27) should be contrasted with con¬ 
ventional thermoelectrics,where the optimal efhciency 
satishes an analogous expression. In conventional ther¬ 
moelectrics, the efhciency of converting heat into electri¬ 
cal energy is limited by the Carnot efhciency rjc- How¬ 
ever, the maximal efhciency which can be reached given 
a set of linear response coefficients is lower than the 
Carnot efficiency by a factor involving the hgure of merit 
ZT = Lil^/deiL, see Eq. (1). In contrast, Eq. (26) de¬ 
scribes the efficiency of converting electrical energy into 
other (e.g., mechanical) forms of energy. This process is 
not fundamentally limited and hence Eq. (26) does not 
contain an analog of the Carnot efficiency. However, it 
still contains an analog of the factor involving the ZT 
value, which contains an appropriate hgure of merit C,. 
Thus, the motor efficiency 77 “°* is bounded by unity, and 
reaches this limit when C —)■ 00 , le., when one of the 
dissipative coefficients Ln or L 33 approaches zero. The 
different form of the hgure of merit, i.e., the absence of 
the coefficients L 13 and L 31 from the denominator, re- 
hects the fact that unlike L 12 and L 2 I) these coefficients 
do not affect entropy production. 


Inserting this into the dehnition (28) of the efficiencies, 
we obtain 


^he ^ ^ Tl-Tr ^ 

Tyhp Tr • 


For heat engines, heat hows from the hot to the cold 
reservoir so that Qr < 0 , while for heat pumps, heat 
hows in the opposite direction, Qr > 0. Thus, we 

hnd that the second law ,5 > 0 implies that the ef- 
hciencies are smaller than the familiar Carnot efficien¬ 
cies, i.e., r]c = {Tr — Tr)/Tr for the heat engine and 
rjc = Tr/{Tr — Tr) for the heat pump. 

The efficiencies for heat engine and heat pumps can 
alternatively be expressed as 


’ Tyhp X 2 J 2 ’ 


(31) 


where the huxes can be expressed through their linear- 
response expressions (16). Maximizing the efficiency as 
for motors and generators, we again hnd Eqs. (24, 25), 
but with the affinity X 2 taking the place of Xi. This 
leads to a maximal efficiency of 


B. Heat engine and heat pump 


Analogous results are obtained when the device is 
driven by a temperature gradient ST at constant chem¬ 
ical potential {Xi = 0), see Fig. 1(b). When the device 
operates as a heat engine, i.e., for L 32 ST/T‘^ < 0 , heat 
hows to the cold reservoir and the system performs work 
on the ac potentials. Conversely, the device operates as 
a heat pump when L 32 ST/T'^ > 0 , where heat is pumped 
to the hot reservoir by the ac potentials. As a result of 
the Onsager symmetry, we have L 23 = —A 32 for time- 
reversal symmetric systems, and we again hnd that the 
second law imposes L 22 > 0 and L 33 > 0 . 

An appropriate measure of the efficiency of a heat en¬ 
gine ? 7 *'® is the ratio of the work per unit time performed 
by the electrons on the ac forces, —W, and the heat leav¬ 
ing the hot reservoir —Qr. (We assume that the left 
reservoir with temperature Tr = Tr — ST is the hot reser¬ 
voir. The efficiency of a heat pump is characterized 
by the inverse ratio. Thus, we have 




he 


1 

T^hp 


W 

Ql 


(28) 


= Vc 


1 + C-i 

1-bC-bi 


with the figure of merit 


c = 


—Lo^^L: 


23-^32 


LooL 


22-^33 


(32) 


(33) 


Equation (32) holds for both heat engines and heat 
pumps, when the appropriate Carnot efficiency rjc is used. 


V. EXAMPLE 

To illustrate these concepts, we consider a quantum 
dot with a single level coupled to two reservoirs with 
chemical potentials /Xq, and temperatures Tq,, a = L,R 
as sketched in Fig. 2. We assume that the dot level and 
the barriers can be modulated periodically in time by ac 
gate potentials. This model can describe a single-electron 
source, similar to the GHz pump realized experimentally 
in Ref. 10. For noninteracting electrons, the model is 
described by the Hamiltonian 

nit) = n,{t) +Tires + nr- (34) 


for the efficiencies of heat engine and heat pump. 

We first show that the second law implies that these 
efficiencies are bounded by the corresponding Carnot ef- 


The first term describes the central conductor which is 
modeled as a discrete chain of N sites with local ener¬ 
gies Em, nearest-neighbor hopping w, and an ac potential 
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Periodic driving 

FIG. 2. Sketch of the device. A single-level quantum dot 
(m = 2) is defined by two tunnel barriers (m = 1,3). They are 
driven by periodic gate potentials Vm{t) = V° cos{ujt -|- Sm), 
with Vi = V§ = 4, V 2 = 23, = 0, ^2 = 7 r/ 2 , and ^3 = 

TT. The tunneling amplitudes between barriers and dot are 
ui = 1 and wl = wr = 0.7 between barriers and reservoirs. 
The barriers and the dot are modeled by a discrete chain of 
A = 3 sites with local energies ei = £3 = 3.3 and £2 = —1 
respectively. The reservoirs have fiR = fi — Sfj, and 

temperature T. 


applied to each site, 


N 




N-l 




+h.c. 


m—1 


m—1 


(35) 

Specifically, we consider a setup with TV = 3 sites, mod¬ 
eling the tunneling barriers (m = 1,3) and the quantum 
dot (to = 2). The site energies (m = 1,2,3) are mod¬ 
ulated by three time-dependent gate voltages of the form 
Vm{t) = cos{ujt + 6m)■ The reservoirs are represented 
by free-electron Hamiltonians for free electrons. 


'^res — E (36) 

oi=L,R,ka 


and tunneling between reservoirs and central system is 
described by 

= - ^ [wadi^Ck^ + h.c], (37) 

(y,ka ,n 


where tIq, denotes the site of the central conductor which 

is in contact with the reservoir a. _ 

The mean charge current Na and heat current en¬ 
tering the reservoir a, as well as the mean power W de¬ 
veloped by the ac forces are calculated within a Floquet 
Green function formalism following Ref. 12, as reviewed 
in App. B. To derive the response coefficients Ly, we ex¬ 
pand the currents Ji = Nr, J 2 = Qr, and J 3 = W/{huj) 
to linear order in huj, see App. C. Explicit expressions for 
the coefficients Lij - in terms of Green functions^^’^®’^® 
or scattering matrices^^ and valid for noninteracting sys¬ 
tems - can be found in App. D. These coefficients can also 
be calculated using an alternative procedure which does 
not rely on the Floquet decomposition, see Refs. 27 and 
28. However, we prefer to use the Floquet approach be¬ 
cause this representation stresses that huj appears in the 




FIG. 3. Maximum efficiency and transport coefficients 
at T = 0 for the motor (M) or generator (G) modes. Inset: 
^max g 0 ];j 0 ]-a,tor/motor with r = —0.7 (/r = 7.2). 


Fermi functions which enter the integrals for the currents 
on the same footing as the chemical potential fj,. This 
provides an alternative argument for identifying hui/T as 
an affinity. 

For illustration, we consider an applied bias Sfi at 
T = 0, i.e., the motor/generator regime, and calculate 
the coefficients listed in App. D for this case. In Fig. 
3, we plot the transport coefficients and the maximum 
efficiency as functions of the chemical potential /r 
of the left reservoir. Large values of the figure of merit 
require a large charge pumping coefficient L 13 along with 
a small value of L 33 L 11 , i.e., low friction or conductance. 
In the absence of driving at the central dot [V 2 {t) = 0], 
the conductance peaks near Ln = I when fj, is in res¬ 
onance with the dot level. Driving the dot level with a 
phase lag relative to the barrier oscillations (^2 — <5^ ^ 0 
for TO = 1,3) favors charge pumping and decreases the 
conductance by dynamically tuning the dot off resonance. 
In this way, high efficiencies can be achieved despite large 
values of L 33 . 

As the chemical potential passes the dot level, 
the pumping coefficient changes sign, and the system 
switches from motor mode [L31611/T < 0; see region M in 
the Fig. 3] to generator mode [L^iSii/T > 0, see region G 
in the Fig. 3]. The efficiency becomes minimal when the 
chemical potential is resonant with the dot level, where 
the conductance is maximal and pumping vanishes by 
particle-hole symmetry. 

The device can also operate as a heat engine or pump 
when imposing a temperature gradient. As this requires 
finite T, quantum effects are less pronounced and effi¬ 
ciencies are lower than those shown in Fig. 3. However, 
we find that for appropriate parameters these may still 































be as high as « OArjc. 


VI. SUMMARY 


Motivated in part by Jarzynski’s equality^^ and 
Crook’s theorem,there has been much interest in quan¬ 
tum thermodynamics, including fluctuation relations, 
work fluctuations, and the thermodynamic description of 
strongly coupled systems.Here, we provided a gener¬ 
alized thermoelectric framework to analyze the thermo¬ 
dynamics of ac-driven nanoscale systems which explic¬ 
itly accounts for the effects of quantum pumping and the 
related nonconservative forces. We identified the addi¬ 
tional flux and affinity through which these forces enter 
the theory and defined generalized Onsager relations for 
the associated response coefficients. This framework al¬ 
lowed us to define appropriate efficiencies and figures of 
merit which describe quantum motors, generators, heat 
engines, and heat pumps. We illustrated our results for 
a simple quantum-pump device. 
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Appendix A: Response coefficients and 
microreversability 


pression as 


Af = Re 


+ 00 /> + 00 7 

/ CLUJ —iujT^ \ 

d^xt Hw) 


= Re 


J —oo J —oo 

/ +00 

-oo 


dud^xi'^' 


Im 


= lim ■ 

o;—>-0 


J,Fj / \ 

Xt (w) 


w 


(A2) 


where we have used that the spectral function 


Im 


■7, Fj / \ 

Xt (^) 


is odd in w, hence Im 




= 0 . 


25 


Analogously, the matrix elements of A(^ can be written 


as 


Af/ = -i 


= Re 




+ 00 





(A3) 


F- F- F F- 

where Xt" ^ i® Fourier transform of Xt" ^ (''”) = 

-ie{T){ Fi(r),i/(0) )f 

The calculation of the conductivity follows the usual 
procedure of the Kubo formula presented in text books. 
We start by considering an extra perturbation due to the 
coupling to an electric field E{t) = dtA{t). In the Fourier 
domain the extra perturbation is H'{u:) = J • E{u})/(iuj), 
which leads to the definition of the dc conductance 


Im 


A“ = lim ■ 

o;—>-0 


J,J f \ 
Xt (w) 


(A4) 


where Xt'"^ i'^) i® the Fourier transform of x/‘^(t) = 
-f6»(r)( J(r),i(0) )t. 

Similarly, evaluating the forces in linear response with 
respect to Sfj. leads to 


The matrix elements entering A/ read 


/ OO 

dt'(t - t')6{t - t'){ F(t),Fj(t') )t 

-OO ^ 

/ oo 

dTTxt^^ir), (AI) 

-OO 


where we have defined the retarded susceptibility 


Xt 


= -i6{T){ 


J^(t),F)(0) )t- In the latter step 


we have stressed that for evolutions with the operator 
U{t) = being Ht = H (V(f)) the frozen Hamil¬ 

tonian, the actual time argument of the integrand of (AI) 
is T = t — t'. Representing the susceptibility in terms of 
the Fourier transform we can also write the previous ex- 


A{" = lim 


Im 


Fj ,J / \ 

Xt (^) 


UJ 


(AS) 


F • J F ■ J 

where Xt^’ i® f^e Fourier transform of Xt^' ('^) = 
-*d(r)([F)(r),J(0)])t. 

The above definitions indicate that the susceptibili¬ 
ties xf"^A with Oi a generic operator, satisfy microre¬ 
versibility with respect to r. It can be directly verified 
that 


Oi,Oi 


Xt 


(-r) = -i6»(-T)(^di(-r),dj(0)J)t = 

= i6»(-r)( 6j(r),6i(0) )t = 

= -i9{-T)f —Im[x/®’°'(a;)]e- 

J — OO ^ 


(A6) 
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Hence, under a transformation r —> —r the coefficient 
Kij transforms to 


Re 


^+oo 


f*+oo 


— dTT9{-T)e- 


Im 


= lim 

Cl »—>-0 


Oj ,Oi / \ 

Xt (w) 


= A' 


w 


o,,o. 


(A7) 


In the last step we have used drre = ;j 2 + 
iTrS'(uj). 

In the presence of a magnetic held B, a time-reversal 
transformation implies changing B ^ —B in the Hamil¬ 
tonian 'Ht dehning the frozen density matrix pt used 
to evaluate the expectation values. This property 
leads to the following Onsager relations for the usual 
susceptibilities in the presence of B, {B,uj) = 

SiSjxf^’^'i—B,ui), where the signs Si, Sj = ± depend on 
the parity of the operators Oi, Oj under a time-reversal 
transformation. 


Appendix B: Green function formalism 


where V(n) are the Fourier components of V(t) = 
being V{t) a matrix with diagonal ele¬ 
ments Vm{t). In the above expressions we introduced the 
hybridization matrix Fq, which has a single element at the 
contact with the reservoir equal to \wa\^2TT^j^^6{E — 
Eka) ■ For practical uses it can be considered in the wide 
band limit, thus, independent of E. The Fermi-Dirac 
distribution fa{E) = [1 -|- characterizes 

the thermal occupation of the electrons in the reservoirs 
(from now on we set the Boltzmann constant = 1). 

The other possible approach is the Floquet scatter¬ 
ing matrix formalism used in Ref. 11. The elements 
Sij{Em, En) of the Floquet scattering matrix s{E), with 
En = E+nhui, are the amplitudes for an electron to scat¬ 
ter from lead j to lead i after acquiring m — n Floquet 
quanta hoj. The general relation between the Floquet 
scattering matrix elements and the Fourier coefficients 
for the Green’s function is the generalized Fisher-Lee 
formula^*^ 

Sij A^) = SijSmn jQij E^')- (H7) 


In Ref. 12 it was shown that the averaged charge Na 
and heat currents entering the reservoir a, as well 
as the mean power W developed by the ac forces can be 
written in terms of the retarded Green function of the 
central structure connected to the reservoirs expanded in 
the Floquet-Fourier transform as 


OO 

n——oo 







(Bl) 


Appendix C: Linear response 


In order to calculate the currents Ji, Z = 1,2,3 up to 
linear order in hui, 6p and ST we perform the following 
expansion of the Fermi function entering the integrands 
of Eqs. (B2), (B4) and Eq. (B6) 


/a(A -f nfuj) - fa{E) + ntiwdEfa{E) 

df{E)^ , df{E){E-p) 
dE dE T 




(Cl) 

T). 


This function is calculated from the Hamiltonian 'H{t) by 
solving the Dyson equation (see Refs. 12 and 40). 

The resulting expression for the charge current is 

^“ = T E [•C/5 in, A), (B2) 


n,/3 


with 


Tapin^E) = \Gapin,E)\'^rpta. 
The heat current reads 


with 


To 

Qa — ^ ol l^a. 5 

e 


A 


(B3) 

(B4) 


A„ = /dA E 0 - [fp iE)-U (A + ntnoWcp {n,E). (B5) 

Similarly, the work performed by the ac potentials can 
be written as 

_ 1 f+o° 

1F = --E / dEntkjfa{E) 

XImjTr \v{n)g{n + I, E)vj^ {I, A)]} , (B6) 


We also evaluate Q{n,E) up to linear order in w by 
expanding the Dyson equation in powers of w (see 35 
and 36) . Up to the Hrst order in w it reads 


G{t,E) ~ G/(t,A) -hiG(i)(t,A), 


(G2) 


with G{t, A) = J2^=-oo E)e The first term is 

the frozen Green function 


Gf{t,E) = 


-1 


lE-Hl + i- 
2 


(G3) 


corresponding to the frozen Hamiltonian at time t, = 
TLcit) (F collects the hybridization functions of the reser¬ 
voirs). The next term is first order in w. It reads 

GW{t,E) = ^dEdtGfit,E)+Ait,E), (G4) 

where 

A=^(dEGfit, A)^G/(t, A)-G/(t, E)^dEGfit, A)J . 

(G5) 
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The expansion of the Floquet scattering matrix up to the 
first order in the driving frequency w cast 


^ij (F/, En) 



nhu) 


2 




dESij {t, E) + hwAij {t, E)]. 


(C6) 


and of the frozen scattering matrix they satisfy 


11 


^ ^ ^ 2 [ dt dE dE dt 


(C7) 


Equation (B7) defines an explicit relation between A and 


Here Sy (t, E) is the frozen scattering matrix. The matrix 
elements Aij{t,E) define a first order correction to the 
adiabatic scattering matrix. The frozen scattering matrix 
Sij{t,E) as well as Aij{t,E) do not change significantly 
on the energy scale hu> and T and depend on the specific 
realization of the scatterer. Anyway, it can be shown 
that, due to the unitarity of the Floquet scattering matrix 


Appendix D: Transport coefficients 

Substituting the expansions for the Fermi function, Eq. 
(Cl), and for the Green function, Eq. (C2), into Eqs. 
(B2, B4) and Eq. (B6) and collecting terms up to first 
order in the affinities Ai = ^, A 2 = p- and X 3 = ^ 
we obtain: 


Lit — — 


T 

hf 


rT 


^ + 00 


/ 

/o J -00 

rp pT 


dE^\Gi^{t,E)\^rLrn 


r-\-oc 


Li2 — L 21 — — ttp / dt 
hT Jo 


Lio — —F.'^i — — 


31 


T 


T 

2 TTh 


rT 


JO 
+ 00 


—00 

T ^+00 
dt / 

J —00 


dE{E-^i)^\Gf^^{t,E)\^tEtR 


E 22 — I dt 


dE^^lm ■ 
dE 


dE(E - ^\G>,^V,E)rr^Ta 




RRd 


hr. 


L23 — —L32 — — 


' —00 


T 

2 'Kh 


fT E + OO 

dt / 


dE (E - fi) Aijn . 


Gf(,.E)E-E^f 


RRy 


L33 = -^ rdtrdE^ReU 




dE 


dGf{t,E)^ dGfHt,E) ^ 


dt 


dt 


Within the scattering matrix formalism the coefficients read 


Lit — — 


T 

W 


T 

Lt 2 = F 2 I = — 


dtj dE^\sRLit,E)\^ 

dE{E - l) f^\sRL{t,E)f 


Li3 — —L,2i — — 


T 


'31 


27r/i 


0 ^ —00 

T ^+00 
dt / 


T - ^ 

^22 — —J— 

hr 


/o 

j\t dE{E-f,)^-^\§RL{t,E)\^^ 




s{t,E) 


dsAt,E) 

di 


J RR 


L 23 — —L 32 — — 
TT 


L33 — — 


SiT^h. 


T 

2 'Kh 

T 

dt 


rdtrdEiE-,)gi^ 

Iq J-oo 
/• + 00 


'0 


dE^Tr 
dE 


dE 

ds{t,E) dsAt,E) 


s{t,E) 


dsAt^E) 

dt 


J RR 


dt 


dt 


(Dl) 


(D2) 


The matrices A in the Green-function language, and A tribute to the coefficient L33. In particular, they appear 
in the scattering matrix version, in principle seem to con- 
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in an integrand of the form 


2 Ini I Ajj (t, £;) 

ij ^ 


ds*,{t.E)\ 

dt j ■ 


(D3) 


However, as shown in Ref. 27, due to the unitary con¬ 
dition of the frozen scattering matrix ss^ = 1 and the 
property (C7) such term vanishes. In fact, it can be also 
written as 


{Tr 

1 = —iTr 

dts^ A — A^dts 

iTr 

J- dts^s 

= 


= ^Tr - S'^dEsdtS^s) (9ts'''s] = 0. 

Zcj 


(D4) 
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